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ABSTRACT 



We study the pure-gauge QCD phase transition at finite temperatures in the 
dual Ginzburg-Landau theory, an effective theory of QCD based on the dual Higgs 
mechanism. We formulate the effective potential at various temperatures by in- 
troducing the quadratic source term, which is a new useful method to obtain the 
effective potential in the negative-curvature region. Thermal effects reduce the 
QCD-monopole condensate and bring a first-order deconfinement phase transition. 
We find a large reduction of the self-interaction among QCD-monopoles and the 
glueball masses near the critical temperature by considering the temperature de- 
pendence of the self-interaction. We also calculate the string tension at finite 
temperatures. 
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It is believed that the hadron physics is governed by the quantum chromody- 
namics (QCD). Although the high energy phenomena as deep inelastic scattering 
are described by perturbative QCD owing to the asymptotic freedom [1], the the- 
ory becomes highly non-perturbative at low energy. The lattice QCD theory [2] 
was then developed for low-energy phenomena and demonstrated again the good- 
ness of QCD. This gigantic numerical simulation method, however, does not tell us 
yet how these low-energy phenomena as quark confinement and chiral symmetry 
breaking take place. We need therefore some effective theory, which incorporates 
the essence of the low-energy QCD physics and at the same time reproduces the 
observables with a few parameters in the effective theory, as the Ginzburg-Landau 
theory in the superconductivity. 

In recent years, the Kanazawa group [3] proposed the dual Ginzburg-Landau 
theory (DGL) as an attractive effective theory of nonperturbative QCD. The DGL 
theory incorporates the QCD monopole fields as essential ingredients for confine- 
ment of colored particles (quarks and gluons). The QCD monopole field has its 
clear origin through the abelian gauge fixing in the non-abelian gauge theory a la 
't Hooft [4]. In this theory, the QCD vacuum is characterized by QCD monopole 
condensation, which provides a mass to the dual gauge field through the dual Higgs 
mechanism and hence the color electric field cannot freely develop in the QCD vac- 
uum. Therefore, the color electric field originated from one colored object should 
be confined in a small vortex-like tube to end at the other colored object [3,5]. 
This corresponds precisely to the Meissner effect of superconductivity. But, here 
the role of the magnetic and the electric fields are reversed and this phenomenon 
is called as the dual Meissner effect. Such a picture for the color confinement has 
been extensively investigated by recent studies [6,7] based on the lattice gauge the- 
ory, and many evidences of QCD-monopole condensation in the nonperturbative 
QCD vacuum have been reported [6]. 

Suganuma, Sasaki and Toki (SST) studied the DGL theory in further detail [5]. 
They formulated the confinement potential of heavy quarks in a natural way and 
obtained the linear potential in the same form as the vortex-like solution in the 
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superconductivity. SST then showed QCD monopole condensation also induced 
the dynamical chiral-symmetry breaking (D%SB) [5,8,9]. Thus, the DGL theory 
describes both the confinement and D%SB of QCD in the non-perturbative region. 
These two nonperturbative features, the confinement and D^SB, would be changed 
in the high-temperature system, which is expected to be realized as the quark- 
gluon-plasma (QGP) in the ultra-relativistic heavy-ion collisions or in the early 
universe [11]. Nowadays, the finite-temperature QCD including the QGP physics 
is one of the most interesting subjects in the intermediate-energy physics [11]. 

In this paper, we would like to develop the thermodynamics of the DGL theory 
[12] and study the change of the properties in the QCD vacuum with tempera- 
tures especially in terms of the deconfinement phase transition. To this end, we 
concentrate on the pure-gauge QCD case, where glueballs appear as the physical 
excitation. Although such a pure gauge system is different from the real world, it 
is regarded as a proto-type of the real QCD and is well studied by using the lat- 
tice QCD simulation [2]. It is worth mentioning that our framework based on the 
DGL theory can be extended to include the dynamical quarks straightforwardly 
[3,5] keeping the chiral symmetry of the system, which is explicitly broken in the 
color-dielectric model [10] or in the lattice QCD with the Wilson fermion [2]. 

The DGL Lagrangian [3,5,12] relevant for the QCD vacuum is written by the 
dual gauge field = (B^, B^) and the QCD-monopole field = 1, 2, 3), 

1 3 
£dgl = —^B v - dvBrf + - 9& ■ B^)Xa\ 2 - A(|xa| 2 - v 2 ) 2 }. (1) 

a=l 

Here, e a denotes the magnetic charge of the QCD-monopole field Xa- e*i = (1, 0), e*2 =| 
(— |, — ^) and e*3 = (— ^, ^). The dual gauge coupling constant g satisfies the 
Dirac condition, eg = 4ir, with e being the gauge coupling constant. The strength 
A for the self-interaction of the QCD-monopole field and the vacuum expectation 
value v are the parameters of the DGL theory. In principle, the values of these 
parameters can be extracted from the lattice QCD data, but it is practically dif- 
ficult at present. Hence, these parameters are determined by fitting to various 



3 



low-energy observables. The DGL Lagrangian (1) is obtained by integrating over 
in the original DGL partition functional in the pure gauge case [3,5]. 

We investigate the effective potential in the path integral formalism. The 
partition functional is written as 




where we take the quadratic source term [13] instead of the standard linear source 
term [1,11]. As is well-known in the 4 theory [1,11], the use of the linear source 
term leads to an imaginary mass of the scalar field Xa in the negative-curvature re- 
gion of the classical potential, and therefore the effective action cannot be obtained 
there due to the appearance of "tachyons". In this respect, there is an extremely 
advanced point in the use of the quadratic source term [13], because the mass of 
the scalar field Xa is always real even in the negative-curvature region of the clas- 
sical potential owing to the contribution of the source J to the scalar mass. [See 
Eq.(6).] Then, one obtains the effective action for the whole region of the order 
parameter without any difficulty of the imaginary-mass problem. Moreover, the 
effective action with the quadratic source can be formulated keeping the symmetry 
of the classical potential. Since this method with the quadratic source term is quite 
general, it is convenient to formulate the non-convex effective potential in the 4> 4 
theory, the linear a model or the Higgs sector in the unified theory [1]. 

The vacuum expectation value of Xa (a=l,2,3) is the same value x due to the 
Weyl symmetry [3] , and therefore we separate the QCD-monopole field Xa into its 
mean field x an d its fluctuation Xa as 

Xa = (x + Xa) exp (if a ). (3) 

Here, the phase variables f a have a constraint, Y2a=i ^ a = 0> where two inde- 
pendent degrees of freedom remain corresponding to the dual gauge symmetry 
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[U(l) 3 x U(l) 8 ] m [3,5]. When QCD monopoles condense, the phase variables £ a 
turn into the longitudinal degrees of the dual gauge field B^, which is the dual 
Higgs mechanism. 

Since we are interested in the translational-invariant system as the QCD vac- 
uum, we consider the ^-independent constant source J, which leads to a homoge- 
neous QCD-monopole condensate. In the unitary gauge, the Lagrangian with the 
source term is rewritten as 



£dgl - J J2 l^l 2 = - 3 Jx 2 - 2x[2A(x 2 - v 2 ) + J]J2x a 

Q=l Q=l 
1 1 3 

- - A {d,B v - d v Brf + -m\Bl + Yj d ^a? - m\xl] (4) 

a—l 

3 

+ Y,i9 2 (^ ■ B^Hxl + nXa) - K^Xxl + xt)}, 

Q=l 

where £ c l(x) is the classical part, 

C cl (x) = -mx 2 -v 2 ) 2 . (5) 

Here, the masses of Xa and B^ are given by 

m 2 x = 2X(3X 2 -v 2 ) + J = 4X X 2 , m 2 B = 3g 2 X 2 , (6) 

where we have used the relation between the mean field X an d the source J, 

J=-2\(x 2 -v 2 ). (7) 

This relation is obtained by the condition that the linear term of Xa vanishes. It 
is remarkable that the scalar mass m x is always real owing to the source J even in 
the negative-curvature region of the classical potential, X < v/y/3. 
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Integrating over B^ and Xa by neglecting the higher order terms of the fluctu- 
ations, we obtain the partition functional, 

Z[J] = exp (i J d 4 x{C cl (x) - 3Jx 2 }) [BetiiD^T^BetiiD- 1 )]- 3 / 2 , (8) 

where the exponents, —1 and —3/2, originate from the numbers of the internal 
degrees of freedom. Here, Db and D x are the propagators of and Xa m the 
QCD-monopole condensed vacuum, 



£J _ | _ MM * JJ _ * /g\ 

m 2 B J k 2 — in'j, -He' x k 2 — m 2 , + ie 



in the momentum representation. Hence, the effective action is given by the Leg- 
endre transformation [1], 

T{X) = -i\nZ[J} + J d 4 x3Jx 2 = J d A xC cl {x) +i lnDet^L^ 1 ) + |i lnDet^D" 1 ). 

(10) 

The functional determinants are easily calculable in the momentum space, and we 
obtain the formal expression of the effective potential [1] , 

V CS (X) = -T(x)/ J d 4 x = 3X( X 2 - v 2 ) 2 + 3 J ln(m| - k 2 - ie) 



i v 11 ) 
3 f d k 2 ,2 ■ \ 

+ 2j J(2^ Hm x~ k ~ l€) 



In the finite-temperature system [11], the partition functional Z is described 
by the Euclidean variables; x$ = —ir, and the upper bound of the r integra- 
tion is /3 — 1/T with T being the temperature. Then, the /co -m tegration in the 
functional determinant becomes the infinite sum over the Matsubara frequency 
[11]. The effective potential at finite temperatures physically corresponds to the 
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thermodynamical potential, and is given by 



V eS (x; T) = 3A(x 2 - t; 2 ) 2 + 3T f; f ln{(2mrT) 2 + A; 2 + m|} 

n=— oo ^ ' 



(12) 



in the DGL theory. Performing the summation over n and the angular integration, 
we obtain the final expression of the effective potential at finite temperatures, 



oo 

V eS (x; T) = 3A(x 2 - ^ 2 ) 2 + 3 J / rfA;A; 2 In (l - e -V^/r) 

o 

oo 

+ 3 -^j dkk 2 In (l-e-v^f/^ 



(13) 







where and m x are functions of \ as shown in Eq.(6). Here, we have dropped 
the T-independent part (quantum fluctuation), because we are interested in the 
thermal contribution to the QCD vacuum [12]. 

We show in Fig.l the effective potential at various temperatures (thermody- 
namical potential), V e s(x] T), as a function of the QCD-monopole condensate x, an 
order parameter for the color confinement. The parameters, A = 25, v — 0.126GeV 
and g = 2.3, are extracted by fitting the static potential in the DGL theory to the 
Cornell potential [5]* These values provide m# = 0.5GeV and m x = 1.26GeV 
at T = 0. The (local-)minimum point of V e s(x;T) corresponds to the physical 
(meta-)stable vacuum state. As the temperature increases, the broken dual gauge 
symmetry tends to be restored, and the QCD-monopole condensate in the physi- 
cal vacuum, Xphys(T), is decreased. A first order phase transition is found at the 
thermodynamical critical temperature, Tq — 0.49 GeV, and the QCD vacuum be- 
comes trivial, Xphys(^) = 0, for T > Tq- This phase transition is regarded as the 



* We examined several possible parameter sets, and found a small parameter dependence on 
our results shown in this paper. 
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deconfinement phase transition, because there is no confining force among colored 
particles in the QCD vacuum with Xphys(^) — 0. 

We show the behavior of the QCD-monopole condensate in the physical vac- 



0.126 GeV at T = 0, and the QCD-monopole condensate decreases monotonously 
up to Xphys — 0.07 GeV at the upper critical temperature T up = 0.51 GeV, where 
the minimum at finite x disappears in V e g(x;T). On the other hand, the local 
minimum is developed at x = in V e g(x;T) above the lower critical temperature 
^low — 0.38 GeV, which is analytically obtained by using the high-temperature 
expansion [11,12], 



The minimum value of V e s(x;T) at x = becomes deeper than that at finite x 
above the thermodynamical critical temperature Tc = 0.49 GeV. Here, we get the 
first-order phase transition because we have considered full orders in x 2 as shown 
in Eq.(13). On the other hand, Monden et al. [12] did not get the first-order phase 
transition due to the use of only the lowest order in x 2 m the high-temperature 
expansion [11], and therefore they had to introduce the cubic term in Xa m the 
Lagrangian. 

Here, we consider the possibility of the temperature dependence on the param- 
eters (X,v) in the DGL theory. The critical temperature, Tq = 0.49 GeV, seems 
much larger than the one of the lattice QCD prediction, Tq — 0.2GeV [2]. How- 
ever, we should remember that the self-interaction term of Xa has been introduced 
phenomenologically in the DGL Lagrangian. In particular, the asymptotic free- 
dom behavior of QCD leads to a possible reduction of the self-interaction among 
QCD monopoles at high temperatures. Hence, we use a simple ansatz for the 
temperature dependence on A, 



uum, Xphys(T), as a function of the temperature T in Fig. 2. One finds Xphys 




(14) 




(15) 
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keeping the other parameter v constant. Here, the constant a is determined as 
a = 0.96 so as to reproduce T c = 0.2GeV. (We take A(T) = for T > T c /a.) The 
results for the monopole condensate Xphys(^) are shown in Fig.3. The qualitative 
behavior is the same as in the above argument with a constant A. We find a 
weak first-order phase transition in this case also. Here, we find a large reduction 
of the self-interaction of the QCD monopoles near the critical temperature Tc: 
A(T ~ Tc) ~ 1 is considerably smaller than A(T = 0) = 25. It would be an 
interesting subject to examine such a large reduction of A(T) near Tc from the 
study of the monopole action in the lattice QCD [7]. 

Next, we investigate the variation of the masses of the dual gauge field and 

— * 

the QCD-monopole field Xa at finite temperatures. Here, and Xa would appear 
as the color-singlet glueball field with 1 + and + , respectively [3,6,8]. The glueball 
masses, and m x , at the finite temperature T are given by 

m B (T) = V3g Xph y S (T), m x (T) = 2v^(T)Xphys(^) (16) 

as shown in Eq.(6). In Fig. 4, We show ms(T) and rn x (T) as functions of the 
temperature T using variable A(T) in Eq.(15). (Their behaviors are almost the 
same as the case of a constant A except for the difference of the value of Tc-) It is 
worth mentioning that mg(T) and m x (T) drop down to m#,m x ~ T^(=0.2GeV) 
from m5,m x ~ 1 GeV near the critical temperature Tc- In other words, the 
QCD phase transition occurs at the temperature satisfying rriB,m x ~ T, which 
seems quite natural because the thermodynamical factor l/{exp(\/k 2 + m 2 /T)±l} 
becomes relevant only for m <>T. Thus, our result predicts a large reduction of the 
glueball masses, tub and m x , near the critical temperature Tc- It is desirable to 
study the change of the glueball masses at finite temperatures, especially near Tc, 
in the lattice QCD simulation with the larger lattice size and the higher accuracy. 

We investigate the string tension k at finite temperatures, since k is one of 
the most important variables for the color confinement, and controls the hadron 
properties through the inter-quark potential. We use the expression of the string 
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tension k(T) provided by SST [5], 



k(T) 



e 2 m 2 B {T) 



24tt 



In 



m 2 B (T) + m 2 x (T) 



(17) 



where the glueball masses mg(T) and m x (T) are given by Eq.(16). The results 
are shown in Fig. 5 as a function of the temperature T. In the case of constant 
A, the string tension k(T) decreases very gradually up to the temperature, T up = 
0.51 GeV. On the other hand, in the case of variable A(T), the string tension k(T) 
decreases rapidly with temperature, and k(T) drops down to zero around Tq = 
0.2 GeV. Hence, one expects a rapid change of the masses and the sizes of the 
quarkonia according to the large reduction of k(T) at high temperatures. We plot 
also the results of the lattice QCD simulation in the pure gauge [15] by black dots 
near and below the critical temperature, Tc = 0.2 GeV. We find our results with 
variable A(T) agree with the lattice QCD data. 

We discuss further the temperature dependence of the parameters (A,i>) in the 
DGL theory. Definitely, we should follow the lattice QCD data for this determi- 
nation as the case of the Ginzburg-Landau theory of superconductors extracting 
the temperature dependence from experiments. Since there exists the lattice QCD 
data on the string tension k [15], we try to reproduce k by taking a simple ansatz 
on A and v. We try the following ansatz, 



where the constant a is determined so as to reproduce Tq = 0.2GeV. (We take 
B(T) = for T > Tc/a.) The variable B(T) corresponds to the bag constant, the 
energy-density difference between the nonperturbative vacuum (\xa\ 0) an d the 
perturbative vacuum (\xa\ — 0) in the DGL theory; see Eq.(13). The ansatz (18) 
suggests the reduction of the bag constant at high temperatures, which provides 
the swelling of hot hadrons by way of the bag-model picture. Since we have already 
examined a typical case for variable A(T) keeping v constant, we show here another 




(18) 
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typical case for variable v(T) keeping A constant. The string tension k(T) in the 
variable v(T) case with a = 0.97 is shown by the dashed line in Fig. 5. We find 
almost an identical result and find again a good agreement with the lattice QCD 
data. Other combinations on A(T) and v(T) under the relation (18) also provide 
equally good results on k(T). 

Finally, we investigate the relation between the scalar glueball mass m x (T) 
and the string tension k(T). For variable A(T) keeping v constant, one finds from 
Eq.(17) an approximate relation, 

mJT) (24tt)V2 , x 

^1— ~ - - — ~ 1.6, 19 

near the critical temperature Tq- On the other hand, for variable v(T) keeping A 
constant, the glueball masses at finite temperatures, mg(r) and m x (T), are shown 
by the dashed line in Fig.4, and Eq.(17) leads to a simple relation, 

mJT) ( 2A \ 1/2 

-3.0, (20) 



JUT) Wln{(32 2 + 4A)/3£ 2 } 

for the whole region of T. Thus, the DGL theory suggests a proportional relation 
between the scalar glueball mass and the square root of the string tension at least 
near Tq. It is worth mentioning that Engels et al. [14] obtained a similar relation, 
m GR (T) = (1.7 ± 0.5)a//c(T), for the lowest scalar glueball at finite temperatures 
from the thermodynamical study on the SU(2) lattice gauge theory. Eqs.(19) and 
(20) can be examined from the thermodynamical study on the glueball mass in the 
lattice QCD, which may also reveal T-dependence on the parameters in the DGL 
theory. 

We have studied the properties of the pure-gauge QCD vacuum at finite tem- 
peratures in the dual Ginzburg-Landau (DGL) theory, where the color confinement 
is realized through the dual Higgs mechanism. We have formulated the effective po- 
tential at finite temperatures (thermodynamical potential) using the path-integral 



11 



formalism. We have used the quadratic source term instead of the linear source 
term. The use of the quadratic source term overcomes the problem of the imagi- 
nary scalar mass, which is encountered in the case of the linear source term as is 
well-known in the 4 theory. 

We have found the reduction of the QCD-monopole condensate at finite tem- 
peratures, and have found a first-order deconfinement phase transition at the crit- 
ical temperature Tq — 0.49GeV using the temperature-independent parameters. 
The QCD-monopole condensate vanishes and the broken dual gauge symmetry is 
restored above Tq- We have considered the temperature dependence of the QCD- 
monopole self-interaction noting Tq = 0.2GeV as the lattice QCD simulation in- 
dicates. We have found a large reduction of the QCD-monopole self-interaction 
near the critical temperature. We have investigated the temperature dependence 
of the glueball masses, tub and m x , and have found their large reduction near the 
critical temperature T<j: rag, m x ~ Tq. We have calculated also the string tension 
at finite temperatures. The results agree with the lattice QCD data both in the 
variable A(T) and in the variable v(T) cases. 

In particular, the glueball mass reduction at high temperatures would be an 
important ingredient in the QCD phase transition. In the pure gauge, there are only 
glueball excitations with the large masses (<: lGeV) at low temperatures [2,14], and 
therefore it seems unnatural that the QCD phase transition takes place at a small 
critical temperature, Tq ~ 0.2GeV. This problem would be explained by the large 
reduction of the glueball mass near the critical temperature as is demonstrated 
in this paper. In other words, this glueball-mass reduction may determine the 
magnitude of the critical temperature Tc in the QCD phase transition. 

It is also very interesting to investigate the thermal effect of the dynami- 
cal quark in the DGL theory, especially in terms of the chiral phase transition 
[2,11]. The chiral-symmetry restoration as well as the deconfinement phase transi- 
tion is expected to happen in the DGL theory, because our previous works [5,8,9] 
showed the close relation between the color confinement and the dynamical chiral- 
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symmetry breaking in the DGL theory. In the presence of dynamical quarks, there 
would be two important ingredients. One is the glueballs, and Xa, character- 
izing the confinement properties of the QCD vacuum. Their mass would decrease 
near Tq according to the deconfinement phase transition. The other is the light 
hadrons as pions, which have a close relation to the spontaneous chiral-symmetry 
breaking [1]. The thermal contribution of light hadrons as pions would lower the 
critical temperature Tq. The lattice QCD simulations showed the close values of 
the critical temperatures between the pure-gauge case (Tq ~ 0.2GeV) and the case 
with dynamical quarks (Tc — 0.15GeV) [2]. These close values of Tq would be 
explained if the glueball-mass reduction at high temperatures plays the dominant 
role in the QCD phase transition even in the presence of dynamical quarks. Such 
a conjecture on the glueball mass reduction would provide an interesting subject 
in the lattice QCD [14]. 

We would like to thank H. Monden and S. Sasaki for stimulating discussions on 
the DGL theory. One of the authors (H.S.) is supported by the Special Researchers' 
Basic Science Program at RIKEN. 
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FIGURE CAPTIONS 



1) The effective potentials at various temperatures as functions of the QCD- 
monopole condensate % The numbers beside each curve are the tempera- 
tures. The absolute minimum points of the effective potentials are shown by 
crosses. 

2) The QCD monopole condensate Xphys(^) at minima of the effective poten- 
tial as a function of the temperature T. The solid curve denotes Xphys(^) 
corresponding to the confinement phase, which is the absolute minimum up 
to Tq =0.49 GeV and becomes a local minimum up to T up =0.51 GeV. A 
minimum appears at x = above T\ ow =0.38 GeV and becomes the absolute 
minimum above Tq =0.49 GeV. The dot-dashed curve denotes the value of 
X at the local maximum. 

3) The QCD monopole condensate Xphys(^) a t minima of the effective potential 
as a function of the temperature in the case of variable A(T) so as to reproduce 
Tc =0.2GeV. The meanings of the curves are the same as in Fig.2. 

4) The masses of the glueballs at various temperatures: ms(T) and m x (T). The 
solid lines denote the case of variable A(T) with a constant v. The dashed 
lines denote the case of variable v (T) with a constant A. A large reduction of 
these masses is found near the critical temperature. The dotted line denotes 
m = T. The phase transition occurs at the temperature satisfying m#, m x ~ 
T. 

5) The string tension k(T) as a function of the temperature T. The solid and 
dashed lines correspond to the variable A(T) case with a constant v and the 
variable v(T) case with a constant A, respectively. The constant (A, v) case 
is also shown by the thin line. The lattice QCD results in the pure gauge in 
Ref.[15] are shown by black dots near the critical temperature. 
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